Searching for Stoponium along with the Higgs boson 
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Abstract 

Stoponium, a bound state of top squark and its antiparticle in a super symmetric model, may 



be found in the ongoing Higgs searches at the LHC. Its WW and ZZ detection ratios relative to 
the Standard Model Higgs boson can be more than unity from WW* threshold to the two Higgs 
threshold. The 77 channel is equally promising. A stoponium mass 135 to 150 GeV is sererely 
constrained by the ATLAS and CMS experiments. 
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Discovery of a Higgs boson is a top priority of LHC experiments, as is the search for 
the supersymmetry which is an attractive candidate for physics beyond the SM. Here we 
make the observation that there is a possibihty of finding supersymmetry in the LHC search 
for the Higgs boson. 

The stop ti, a scalar superpartner of the top quark, is expected to have the hghtest mass 
of all the squarksjlj], and may plausibly be lighter than the top quark. If the mass difference 
between ti and the LSP neutralino A^^i is small and the tree-level decays, ti t Ni and 
ti — 7- 6 Ci (where Ci is the lightest chargino), are not kinematically allowed, the ti becomes 
long-lived and tl ti will form a bound state called stoponium. This expectation is supported 
by the result [2] that the partial width of ti — c A^^i occurring at loop level is negligibly small 
compared with the binding energy of the stoponium, a few GeV. 

The possibility of stoponium discovery at hadron colliders was considered long agojs, 4|. 
Its production via gg-fusion and its decays are quite similar to the heavy quarkonium of the 
fourth generation quarks j^. The production amplitude is proportional to the wave function 
at the origin, and so the S-wave J^^ = 0^^ bound state, denoted here as a, is expected to 
have the largest production cross section. 

Detection of stoponium is complementary to the detection of stop. Actually current LHC 
data already give a strong limit on the light stop massQ| by considering the possible decay 
channels of — )■ biW* and — tNi. For the parameter space where these decay modes 
are kinematically forbidden and stop is long-lived, as required for stoponium existence, this 
constraint does not apply. The stop decays to charm quark and LSP neutralino at loop-level, 
and this decay mode is notoriously difficult to be identified, because of hadronic effects and 
the small phase space. In this parameter region, stoponium detection may be the best way 
to search for supersymmetry. 

In this Letter we refine the calculation of the production and decays of the a appropriate 
to the LHC experiments at 7 TeV (LHC7) and consider the possibility of finding the a in 
the ongoing LHC Higgs searches. Related work can be found in ref. 7|, and more recently 
in ref. 8], at tree level, and the NLO radiative corrections are considered in refs. jol. [lo| . We 
demonstrate that a could be found in the SM Higgs search in the 77 and W*W* channels 
at LHC7. Stoponium can be distinguished from a Higgs boson by differences in decay 
branching fractions. 
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Stoponium Production Cross Section The production cross section of stoponium a in 
hadron colliders is mainly via gg fusion, similarly to the production of a Higgs boson h^. 
The cross sections are proportional to the respective partial decay widths to gg. The pro- 
duction cross section of has been calculated in NNLO 11|, and by using this result (20I we 
can directly estimate the production cross section of a as 

aipp ^ aX) = aipp -> h'X) x ^^^o'^^^g^ ■ (1) 

By using the r{a — )■ gg) partial width given later and r{h^ — )■ gg) of the SM we can predict 
a{pp —7- aX). The result is compared with the SM Higgs production in Fig. [1] 




FIG. 1. The production cross section of o'[pb] from gg fusion (solid blue), cr{gg fT)[pb], compared 
with that of the SM Higgs with the same mass mf^o = m^- (solid red). The overall theoretical 
uncertainties [11] are denoted by dotted lines. 



The production of a exceeds that of the SM Higgs boson of the same mass m/jo = ms- 
for nia < 230 GeV. This is because the a production from gg fusion has an amplitude from 
the 4-point coupling at tree level, while production is governed by the one-loop diagram 
of the top quark. Our prediction of cr{gg — t- a) in FigH] includes the ±25% uncertainty 
associated with the theoretical uncertainty on a{gg — t- h^). 

Stoponium Decay For the a decay channels a — t- AB, we consider AB = gg, 77,^7, W~^W~, 
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ZZ^ bb, tt, and h'^h'^. Their partial widths are given by the formula 

where Ai represent the free titl AB amplitude, p is the momentum of particle A(or B) 
in the CM system, and -R(O) is the radial wave function of stoponium at the origin. The 
total width r|°* is the sum of these partial widths. Here we omit the LSP neutralino channel 
AB = iViiVi, which is a suppressed decay modeM, M. All the relevant formula are given 

npri 

in the previous works|^, 13, |8[, so we briefly explain here our method and the selection of 
parameters. 

In calculations of the amplitude Ai = — t- AB), the mass of the lighter higgs is 
fixed to m/jO = 125 GeV, and a maximal-mixing 6f: = tt/A between the two stops ti and is 
taken. The value of tan/3 = Vu/vd, the ratio of vacuum expectation values of the two Higgs 
doublets in SUSY, is taken to be 10, and the higgs mixing angle a is also fixed from the 
tree-level formula, tan2Q! = [(( m\+ My) / (m^ — My))ta.n20], with the choice niA = 800 GeV, 
for which tana = —0.103 .[21| 

The contribution from the heavier higgs and the heavier stop t2 to the amplitudes 
are neglected under the assumption that they are heavy. Amplitudes of t-,u— channel t 
or b exchanges are neglected except for the ti exchange in the h^h^ channel. Then, all 
the amplitudes are described by the contact 4-point interaction and/or the s-channel 
amplitudes when they contribute. The Feynman diagrams of our analysis are shown in 
Fig.H 

For the gg decay of a we include the radiative corrections at NNLO by using the K-factor 
from Higgs production 22 1. For decay to 77 a special attention is taken by ref. joj since this 



branching fraction is larger than that of the Higgs boson for which a strong cancellation 
between top-quark and W loops occurs. We use the R^^^ = T^a — t- 77)/r(cr — )■ hadrons) at 
NLO in Ref.jol, equating their r(cr — )■ hadrons) to our r(5" — ?■ gg). Multiplying R^^^ by our 
r{a gg), we obtain r(cr — »■ 77) . The off-shell WW*{ZZ*) channels in the low-mass a 



case are treated following ref. 
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The gg partial decay width of a is proportional[5| to {\R{0)\/m^)'^. We use the non- 
relativistic quark model with the Wisconsin potential, where a potential term in the in- 
termediate range is added to Cornell potential, to obtain the value of |-R(0)| The 
binding energy, a few GeV, is much smaller than m^^ and the stoponium mass is well 



(a) 



(b) 



(c) 



FIG. 2. Dominant generic diagrams of a decay. Thin dashed hnes with arrows represent the initial 
titl of a and two wavy hnes represent the final states XX. Diagram (a) is taken into account in 
XX = gg,jj,Z'y, diagram (b) in XX = bb, and diagrams (a) and (b) in XX = WW,ZZ, while 
all three diagrams contribute to XX = hPhP. 
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FIG. 3. Decay Branching Fractions of a versus m5-(GeV) for k = 3 and -1. rrij^o is taken to be 125 
GeV. 



approximated by 2mf^. 



With these simplifications, the results depend only upon two quantities, rria and the 
hHit\ coupling A^ojj., denoted as mi^\2 in our previous work^. The is determined 



by a dimensionless parameter k . 23| In our illustrations we consider two values of k, k = 3 
and -1, which correspond to the strong and weak h^tit\ coupling, respectively. The is 
taken as a free parameter varying a wide range, 135 GeV < < 300 GeV. Below 135 
GeV, the mixing between a and (with m/jO taken to be 125 GeV here) may be important, 
although this mixing can be handled if necessary. 

The decay branching fractions of a are given in Fig. [31 In the case of k = 3 [larger 
^H^tii.^ 396) GeV case], the branching fractions of WW, ZZ,h'^h'^,bb,ti are larger than 
those for K = — 1 [smaller Xh"iii*{= 169) GeV case] because of larger contributions from 
the s-channel diagram. WW, ZZ, h^h^ branching fractions become even larger in the 
K = 10 case, where WW is dominant in 2m\Y < < 2m/jO. In the k = — 1 case, the 
decay amplitude to vanishes when ~ 370 GeV because of a destructive interference 
among the 4-point interaction, t-,M-channel ti exchange and the s-channel diagrams. The 
existence of this cancellation has not been previously noted in the literature. 

The total width of a in the mass range = 135 ~ 250 GeV is fairly small, r^°* ~ 2 
to 7 MeV, and in the mass range = 260 ~ 300 GeV, it is at most ~ 40 MeV for the 
K = 3 case. Thus a a resonance would be observed with the width of the experimental 
resolution. At > 2mw the r(a — )■ WW) partial width is negligibly small compared with 
r(/i° — )• WW), and thus a production via vector boson(Vriy, ZZ) fusion is negligible at the 
LHC. 

The ZZ to WW ratio of the a decay branching fraction is predicted to be 0.32~0.36 (for 
K = 3) in the mass range 200 GeV< m^- < 300 GeV, as compared with 0.36~0.44 of the 
SM /i^ in the same mass range. This ratio can be used to check if an observed resonance is 
actually stoponium or not. 

Stoponium Detection compared to SM Higgs Next we consider the detection of a in 
W^W~ , ZZ and 77 channels. The a search can be made in conjunction with the Higgs 
search. The properties of hP at the LHC are well known, so we use them as benchmarks of 
the search for a. 

The a detection ratio {DR) to in the XX channel is defined 15|] by 



tot 



= r r T^tot ' ("^J 

where XX = W^W~, ZZ, and 77. The DR are plotted versus = m/jo in Fig. IHfor the 
two cases k = 3 and -1. 
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FIG. 4. a Detection Ratio (DR) to the SM higgs /i° of Eq. ^ for the XX = W+W~{soM blue), 
ZZ {dashed green), and 77(soHd black) final states for k = 3,-1 versus m5-(GeV). 



In the case of k = 3 (the large A^of^f* case), the a to detection ratio is relatively large in 
both WW and ZZ channels. The ratio is more than 0.5 in the mass range 160 < < 250 
GeV between the WW threshold and the h^h^ threshold. For < 190 GeV the ratio is 
nearly unity. Even in the k = —1 case, the detection ratio is large in the 160 < < 250 
GeV mass range. 

The cross-section of a putative Higgs-boson signal, relative to the Standard Model cross 
section, as a function of the assumed Higgs boson mass, is widely used by the experimental 
groups to determine the allowed and excluded regions of m/jO. By use of the DR in Fig. HI 
we can determine the allowed region of from the present LHC data. Figure [5] shows the 
95% confidence level upper limits on Higgs- like a signals decaying into XX versus for 
XX = WW and ZZ combined (ATLAS [16| [a^) and 77 (ATLAS [ij]). 

For K = 3, is excluded by ATLAS data over wide ranges of 155-227 GeV, while for 
the K = — 1 some regions of are excluded. Similar results are found from CMS datafl?!. 

The a search is also applicable to the Tevatron data. The CDF and DO experiments 
excluded the SM Higgs with mass 158 GeV < niho < 175 GeV from the data of WW, ZZ 
channels. The same data excludes a in the k = 3 case in the mass range, 160 GeV < < 




FIG. 5. the 95% confidence level upper limits of (l/DR) x {acxp/(^igg ^ ^X)). This is the 



signal of a boson decaying into XX relative to the stoponium cross section Wigg — ?> o" — ?> XX) 



18|) data. The cases 



a{gg ^ ^ XX) x DR] for XX = VFVF (ATLAS [16|) and 77 (ATLAS 
K = 3(solid blue) and -1 (solid green) are shown. Similar results for the SM Higgs boson are also 
given(red thin-solid curve). 



177 GeV. 

The 77 final state is very promising for a detection, because the a to detection ratio 
is generally very large in all the mass range of m^-, as shown in Fig. HJ From the 77 data of 
ATLAS the region of stoponium mass 135 < < 150 GeV is already excluded in both the 
cases of k. The ^-dependence of DR in 77 is small below WW threshold because the partial 
widths of all the allowed two-body final states are independent of k. For > 150 GeV, 
the 77 signal of hP is too small to be detected, but the data in this region can determine 
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the existence of a. 

Concluding Remarks We have investigated the possibihty of finding stoponium a at LHC7. 
In the optimistic case of the stoponium mass and couphng, a will be discovered in the 
WW, ZZ, and 77 channels in the search for the SM Higgs h^. The detection rates can be 
comparable to that of the SM h^, in the mass region m^- ~ 160 GeV up to 2m/jO as shown 
in Fig. m The 77 search channel is particularly promising since a detection relative to 
is very large (more than 3) in all mass regions. A stoponium mass in the 135-150 GeV is 
already excluded in a wide range of supersymmetry parameters from the present ATLAS 77 
data. The h^h^ decay of a is another possible mode for discovery, especially for large k. 
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